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Abstract 

We calculate the shear viscosity (77) to entropy density (s) ratio 77/s of a gluon plasma in kinetic 
theory including the gg — > gg and gg — > ggg processes. Due to the suppressed contribution to 77 
in the gg — ► gg forward scattering, it is known that the gluon bremsstrahlung gg <-> ggg process 
also contributes at the same order (0(aj 2 )) in perturbative QCD. Using the Gunion-Bertsch 
formula for the gg — > ggg matrix element which is valid for the limit of soft bremsstrahlung, we 
find that the result is sensitive to whether the same limit is taken for the phase space. Using the 
exact phase space, the gg <-> ggg contribution becomes more important to 77 than gg — > gg for 
a s > 2 x 10~ 3 . Therefore, at a s = 0.1, 77/s ~ 1.0, between 2.7 obtained by Arnold, Moore and 
Yaffe (AMY) and 0.5 obtained by Xu and Greiner. If the soft bremsstrahlung limit is imposed 
on the phase space such that the recoil effect from the bremsstrahlung gluon is neglected, then 
the correction from the gg «-> ggg process is about 10-30% of the total which is close to AMY's 
prediction. This shows that the soft bremsstrahlung approximation is not as good as previously 
expected. 



USTC-ICTS-09-10 



I. INTRODUCTION 



One of the most surprising discoveries at the Relativistic Heavy Ion Collider (RHIC) is 
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that the hot and dense matter (believed to be a quark gluon plasmaJQG 
for reviews) formed in collisions appears to be a near-perfect fluid 
131 ] . The remanent of the non-central collisions shows collective motion (elliptic flow) with 
a shear viscosity (77) to entropy density (s) ratio 77/5 = . 1 ± . 1 (theory) ±0.08 (experiment ) 



[l4 ] . This 77/s ratio is close to a conjectured minimum bound 1/47T Il5lj , which is motivated 

u, y. 



by uncertainty principle 



161 ] and gauge/string duality 
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If 



rj/s implies stronger particle interactions, contrary to the conventional picture, the QGP 
produced at RHIC tends to be a strongly interacting fluid instead of a weakly interacting 
gas. 

However, a recent perturbative QCD calculation of 77/s of a gluon plasma by Xu and 
Greiner (XG) 21| indicates that the gluon elastic scattering gg — > gg does not give 
the dominant contribution. They found that rj/s for the gluon bremsstrahlung process 
99 999 is about 1/7 of that for gg — > gg, which means the contribution to the shear 
viscosity from gg <-> ggg is 7 times as important as that from gg — > gg. This would bring 
rj/s down to 1/Air when strong coupling constant a s ~ 0.6. This implies that the near- 
perfect QGP might still be described by perturbative QCD and that the conventional 
picture could still be valid. Their conclusion is quite different from an earlier study by 
Arnold, Moore and Yaffe (AMY) 22J (for a recent review, see, e.g., 23j). AMY found 
that gg <-> ggg only contributes at 10% level for the three flavor quark diffusion constant 
for a s < 0.3. For comparison, XG have 77/s ~ 0.5 at a s = 0.1, while AMY have rj/s ~ 2.7 
(note that only r\ was computed in 22j], the free gluon s is inserted by us for comparison). 

Both approaches of XG and AMY are based on kinetic theory. However, the main 
points of differences are: 1) A parton cascade model 24J is used by XG to solve the 
Boltzmann equation. Since the bosonic nature of gluons is hard to implement in real 
time simulations in this model, gluons are treated as a Boltzmann gas (i.e. a classical 
gas). For AMY, the Boltzmann equation is solved in a variation method without taking 
the Boltzmann gas approximation. 2) The Ng <-> (N + l)g processes, N = 2, 3, 4 . . ., are 
approximated by the effective g <-> gg splitting in AMY where the two gluons are nearly 



collinear with a small splitting angle, while the gg <-> ggg process is used in XG where 
the bremsstrahlung gluon is soft but it can have a large splittin g an gle with its mother 



for the gg 



999 



gluon. More specifically, in XG, the Gunion-Bertsch formula 
matrix element squared in Eq. (112j) is used. This formula is valid for the limit of soft 
but not necessarily collinear gluon bremsstrahlung. For the phase space, XG uses the 
exact phase space for the three gluon configurations (called "three-body-like" phase space 
in this paper). In principle, if the soft bremsstrahlung limit is a good approximation of 
the gg — > ggg process, one should be able to impose the same limit to the phase space 
as well and get approximately the same result. In this limit, the recoil effect from the 
bremsstrahlung gluon is neglected, and the phase space (for the two hard gluons) is called 
"two-body-like" here. 

In this paper, we will perform a third independent calculation for comparison. We will 
use the same inputs on the Gunion-Bertsch formula for the gluon scattering amplitudes 
(modulo a factor 2 in Eq. (JT2]) ) with the soft gluon bremsstrahlung approximation as XG 
but we will solve the Boltzmann equation variationally as AMY without taking the Boltz- 
mann gas approximation. We will also test the robustness of the soft gluon bremsstrahlung 
approximation by comparing the results with the two- and three-body-like phase space. 



II. KINETIC THEORY 



Using the Kubo formula, t] can be calculated through the linearized response function 
of a thermal equilibrium state 



di' / dt f dx 3 {[T ij (0),T ij (x,t)]), 

J — oo J 



where T iJ is the spatial part of the off-diagonal energy momentum tensor. In a leading 
order (LO) expansion of the coupling constant, there are an infinite number of diagrams 



261 ] . However, it is proven that the summation of the LO diagrams in a weakly coupled 



theory 26( or in hot QED 



271 ] is equivalent to solving the linearized Boltzmann equation 



with temperature-dependent particle masses and scattering amplitudes. The conclusion 
is expected to hold in weakly coupled systems and can as well be used to compute the LO 



transport coefficients in QCD-like theories 



3sH 



and weakly coupled scalar field theories 
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28l| . hadronic gases [29 
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The Boltzmann equation of a hot gluon plasma describes the evolution of the color 



averaged gluon distribution function / = f(x,p,t) = f P (x) (a function of space, time and 



momentum) as 
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Cd*fp = C[f], (2) 



where E p = p for massless gluons. The driving force for the evolution is the particle 
scattering in the microscopic theory described by the collision term C which is a functional 
of /. It is known that, to compute 77 to LO in the coupling constant a s , we need to include 



99 ~ * 99 an d 99 999 scattering in C 
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42]. We will show this more explicitly later. 
In thermal equilibrium, the gluon distribution /° is static, homogeneous and isotropic 
and hence d^f® = at every x. This implies C[f°] = or detailed balance whose solution 
is just the Bose-Einstein distribution function = l/(e Ep ' T — 1). When the system 
is not in thermal equilibrium, there will be momentum flow due to the breakdown of 
detailed balance. The momentum flow can be characterized by a velocity field V(x). The 
deviation from thermal equilibrium can be characterized by the inhomogeneity of V(x) 
or the derivative expansions of V(x). For simplicity, we work in the comoving frame of 
the fluid element at point x with V = and to the order of first derivatives of V. Thus 
the distribution function can be parametrized as 

U = f^-X P {l + f% (3) 

where 

Xp = [A(p)V ■ V + B{p)p {i pfi V^-]] /T , (4) 

and where the symmetric traceless combinations puPj] = p % pP — 5^/3 and VfiV^-i = 
(VjVj + WjVi) /2 — V • VSij/3. Note that the time derivatives do not appear because 
they can be related to the spatial derivatives by virtue of the conservation of energy mo- 
mentum tensor. Analogously the deviation of the energy momentum tensor away from 
its equilibrium value can be parametrized by the bulk (() and shear (rj) viscosities 

STij = C<%V • V - 277V [i ^ ] . (5) 

Using the definition in kinetic theory T^ v = N g J j0p 1 ^ 1 § ]L fp(x)i one obtains 

where N g = 16 is the gluon polarization and color degeneracy. 



Following the standard procedure (see e.g. 26|) and making use of the Boltzmann 
equation satisfied by B(p), Eq.Q can be recast into 



^ /II J^^\Mn^\ 2 ^r6\ Pl +p 2 -p 3 -p,)(l + A°)(l + f 2 )M 



NH fyj d3p * 

80 J f \ (2n)*2Ei 
x [BijiPi) + Bijipz) - Bijfa) - Bijipi)] 2 
N 2 J f" rf 3 p. 



X [B tJ ( P5 ) + By(p 4 ) + fly (p 3 ) - flii(p 2 ) " %(Pl)] 2 , (7) 

where Bij(p) = B(p)(p % p> — and M 12 ^ 34 and M 12 ^ 345 are amplitudes for (7(7 — > 
gg and (7(7 — ► £7 (7^7 processes (or called 22 and 23 processes in this paper), respectively. 
One useful observation is that the right hand sides of Eqs. ([6]) and (GO) correspond to 
integrations over both sides of the Boltzmann equation, or equivalently, a projection 
of the Boltzmann equation. It is certainly easier to solve the projected equation than 
the Boltzmann equation itself. However, there would be an infinite number of solutions 
satisfying the projected equation, even though the true solution is unique, corresponding 
to that which gives the largest 77 (see e.g. 22J). This makes solving 77 a variational 
problem. 

To solve for B(p), we assume it to be a smooth function which can be expanded in a 
set of orthogonal polynomials, 
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B(z) = z y J2 b rB {r) (z)i (8) 

r=0 

where z = P\p\, B^ r '(z) is a polynomial up to z r and the overall factor z y will be chosen 
by trial and error to get the fastest convergence 3jJ. The B^ r \z) polynomials can be 
constructed using the condition 

/ Wf f ° {1 + / P °)IPl^ (r) (^)fl (s) (^) = T% s . (9) 

One can solve the coefficients b r by equating Eqs. §6§ and ([7j). Then, 77 is just proportional 
to b according to Eqs. ([6]) and (Q. For practical reasons, one uses the approximation 
B{z) = z y ^"=0 b r B^'\z) where n is a finite, positive integer. It can be proved that 77 is 
an increasing function of n. Thus, one can systematically approach the true value of 77. 
For y — 1, the series converges rapidly. From n = 2 to 3, rj only changes by ~ 1%. 

In vacuum, |Mi 2 ^ 34 | 2 = (127rct s ) 2 (3 - tu/s 2 - su/t 2 - st/u 2 ) /2 (see e.g. [43j). In 
medium, s = 0(T 2 ). The most singular part of |M 12 ^ 34 | 2 comes from the colinear region 



(i.e. either t ~ or u « 0) which can be regularized by the Hard-Thermal-Loop (HTL) 
dressed propagators for gluons. However XG only used the Debye mass txib = (A-Ka s ) l l 2 T 



as the regulator just as done in Ref. 



441 ] . so for the sake of comparison between AMY 



and XG, we also use mo as the regulator for soft and collinear divergences in this paper. 
We will use the HTL gluon propagators in one of our future study. Thus, we consider the 
near collinear approximation 

(12tt«^ 2 



M- 



12^34 1 



(su/t 2 + St/U 2 )\ ^0 
or «rjO 



(10) 



In the center-of-mass (CM) frame, we can use the crossed symmetry between the u- and 
t-channels and just use two times of the forward angle, t-channel contribution for the sum 
of the forward (t-channel) and backward (it-channel) angle contributions 

,:2 



|M 12 ^ 3 4| 2 « -(12vra s r su/r 



ItfaO 
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where qx is the transverse (with respect to p\) component of q = p 2 — p±. Because small 
qx could also come from large |g 2 | through the w-channel, it is important to note that 
when using Eq. (fTTj) to calculate the collisional integral, we only pick up the near forward 
scattering (around t = q 2 ~ 0) to avoid double counting. 

For the gg — > ggg process, we will take the approximation that the bremsstrahlung 
gluon is very soft (zero rapidity limit) and y/s is much bigger than all transverse momenta. 
Then the exact result of Ref. 43] reduces to the Gunion-Bertsch formula 25] , 



IM- 



12^3451 
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perm(3,4,5) 



(127ra s 



(q T + ml) 2 \4[(q T -k T ) 



m 



(12) 



where we have inserted the regulator m 2 D as in Ref. 44|. Here is the transverse 
component of the bremsstrahlung gluon momentum (p^) and is still the transverse 
component of q = p 2 —p^- The three final state gluons are identical particles. Thus, there 
are 3! permutations of (p3,P4,Ps), each gives the same contribution. As explained above 
in the 22 case, we need to be careful about using the q<r variable. Small qx could mean 
either the forward ( t ~ ) or backward (u ~ 0) scattering. In the convention adopted 
for Eq. (I12p . one can only pick up the near forward scattering (around t = q 2 ~ 0) but 
not the backward scattering otherwise double counting will happen. Our |Mi2^345| 2 is 



derived from the exact result of Ref. 



431 ] . where Lorentz invariant Mandelstam variables 



are used so there is no this ambiguity, after taking the soft bremsstrahlung limit. Eq. 



(I12p is also consistent with the Gunion-Bertsch formula 25| as explicitly demonstrated 
in App. [A] Effectively, the above treatment of collisional inte gral 



s leads to a factor 2 

difference in the gg — > ggg contribution to t] from that of XG 45| (ours is one half of 
XG's). 

Naively the gg — > gg collision rate is oc J dq^\Mi 2 ^^\ 2 = 0(a s ) and the gg — > ggg 
rate is oc J dq\dk\\ Myi^^\ 2 = 0(a 2 ) (as will be discussed below, k\ has an 0(a s ) 
infrared (IR) cut-off). Thus, the gg — > gg process seems more important than gg <-> ggg. 
However, this is incorrect. In gg — » gg, the amplitude is the largest in the forward and 
backward scatterings. But there is no contribution to rj in these cases since there is no 
momentum redistribution. Mathematically, we have the additional suppression factor 
[Bij{pi) + Bij(p 3 ) — Bij(p 2 ) — Bij(pi)} 2 ~ 0(q^) in Eq. (JTj), while no similar suppression 
in gg <-> ggg. Thus, the gg — > gg collision rate is proportional to J dq^M^^ tS 1 q \ 
0(a 2 log a s ), which is of the same order as 0(a 2 ) of gg <-» ggg, up to a logarithm 

This power counting can be used to argue that other processes such as ggg —> ggg and 
99 ~ > 9999 (called 33 and 24 processes) are higher order under the assumption that the 
most important contribution to rj comes from the configurations with at most two hard 
gluons in the initial or the final states. Under this momentum configuration, one observes 
in Eq. (1121) that adding a soft gluon to the 22 process yields a factorizable form for the 
23 matrix element squared. Schematically, 

l -^ 2 ^0(a s p^), (13) 

where pt denotes the small momentum scale with px — O(gr) — 0{kx)- Analogously, 
adding a soft gluon to the 23 process yields 



2n 



1^33(24) 



2 

- ~ 0(a s pr 2 ). (14) 



|M 23 | 2 

Thus, the 33(24) collision rate is smaller than that of 23 by a factor of 
J c?p 2 "|M 33 (24) | 2 / 1 A^23 1 2 = 0(a s \oga s ). This argument can be generalized to other pro- 
cesses as well. Thus, 22 and 23 are the only processes in the LO under this assumption. 

The phase space of the 3-gluon state plays an important role in the collisional integral 
in Eq. ([7]) for gg <-> ggg, which is controlled by the delta-functions for energy-momentum 
conservation. Since we use the Gunion-Bertsch formula, Eq. (I12p . which is valid for soft 
gluon bremsstrahlung, it is consistent to apply the same condition for energy-momentum 
configuration of the 3-gluon state. This can be done by neglecting the recoil effect due to 



the soft gluon bremsstrahlung, i.e. neglecting the momentum of the soft gluon inside the 
delta- functions as is done in App. [A] Therefore, the phase space for the two near collinear 
gluons in 3-gluon state is 2-body-like. Additionally the exact phase space is 3-body-like if 
the momentum of the soft gluon is kept and treated in equal footing as the other gluons 
in the delta-functions. We will see that using the 3-body-like or 2-body-like phase space 
makes a significant difference in the shear viscosity. 



III. LEADING-LOG RESULT 



In the leading-log (LL) approximation, one just needs to focus on the small qr contri- 
bution from the gg — > gg process. After performing the small qx expansion to Eq. ([7]), 
we obtain (g 2 = kna s ) 

Vll -27.1 T " (15) 

T Mi/0) 



which agrees with that of 28j] very well. Using the entropy density for non-interacting 
gluons, s = N g ^-T 3 , we obtain 

& ~ 3 ' 9 , (16) 
s g'Hl/g) 1 } 

This will be used to check our numerical result later. In contrast, we take the Boltzmann 
gas approximation (/° = e~ Ep ^ T ) used by XG, we get t\ll — 44.7T 3 g~ 4 ln _1 (l/g) and 
s = Ng^T 3 , which would yield tjll/s ~ Q.9g~ 4 \n~ 1 (l/g). Thus, the error from taking 
the Boltzmann gas approximation for the LL result of rj/s is ~80%, where ~70% comes 
from 7] and ~10% comes from s. This suggests that the quantum nature of gluons could 
play an important role on transport coefficients, even though they might not be important 
or thermodynamic quantities. In weak coupling regime, e.g. a s = 10~ 3 , the XG result in 



2l| gives r]22/s ~ 5.6 x 10 3 while the LL result gives 7/22/s ~ 2 x 10 4 , which shows a factor 
4 difference. But the difference from the LL result can be narrowed in Israel-Stewart 
theory 46]. 



IV. TREATMENT OF gg <-> ggg 



As mentioned above, both gg — > gg and gg «-> ggg are needed to compute 77 to the 
leading order (0(a~ 2 )). For the treatment of the 23 process, we consider three cases, (a) 

n 



with the 3-body-like phase space for three gluons and with the LPM effect as the cutoff 
for the soft gluon; (b) with the 3-body-like phase space and but with as the regulator 
for the soft gluon; (c) with the 2-body-like phase space and with m D as the regulator for 
the soft gluon. 

In ), the scale of the kx cut-off is set by the Landau- Pomeranchuk-Migdal (LPM) 



effect, as in Refs. 



21| and |22j. Ref. [47j gives an intuitive explanation of the LPM effect: 
for the bremsstrahlung gluon with transverse momentum \kr\, the mother gluon has a 
transverse momentum uncertainty ~ \kr\ and a size uncertainty ~ l/|fcr|. It takes the 
bremsstrahlung gluon the formation time t ~ 1 / 1 /zt | ~ E k /\k T \ 2 to fly far enough from 
the mother gluon to be resolved as a radiation. But if the formation time is longer than 
the mean free path l m f p ~ 0(a~ l ), then the radiation is incomplete and it would be 
resolved as gg — > gg instead of gg —> ggg. Thus, the resolution scale is set by t < l m f p - 
This yields the condition |/ct| 2 > Ek/lmfp ~ 0(a s ) which is confirmed through carefully 



derivations in Ref. 



48]. 



Here the mean free path l m fp is given by the collision rate R ~ l// m / p which sets the 
scale of the LPM effect is computed via the detailed balance rate. After integration, the 
Boltzmann equation of Eq.([2]) can be written as 

— = n (Rgain — Rloss) , (17) 

where we have used n = J j^zfp- Then the collision rate is the detailed balance rate in 
thermal equilibrium, 



p — -nequil. _ T3 i D if? 
K = K gain ~ K 22 + -^23 + -«32, 



where 



^ = §/n-^^l M ^34| 2 (2vr)V(p 1+P2 - P3 -P4) 

x/?/ 2 (l + /3 ) (l + /4°), 

R23 = ^/n^^l M -;34 5 ^ 

xA°/2 Tl+/3°) (l + /4°) (l + / 5 °), 
^32 = ^23- (19) 

Note that our definition of R is the same as that of XG. The phase space for three 
gluons is 3-body-like in R23. And, as mentioned above, only near forward scattering 
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FIG. 1: (color online) R22 and -R23 of Eq. (|19|) shown as functions of a s for BE and Boltzmann 
gas. 

(around t = q 2 ~ 0) is included to avoid double counting (see App. IXj) . which gives an 
additional factor 1/2 compared to XG [21]. We have computed R self-consistently since 
R23 also depends on R. Our R22 and R23, together with the results for the Boltzmann 
gas approximation (/° = e~ Ep ^ T , 1 + /° — > 1), are shown in Fig. 1. Our .R22, which uses 
Bose-Einstein (BE) statistics, is close to the Boltzmann gas result. Our R23, however, gets 
an enhancement for a s < 0.04 from the enhancement factor (1 + /°) which is inversely 
proportional to the soft gluon's bremsstrahlung energy. This enhancement in R makes 
the gg <-> ggg contribution to rj smaller in the BE case than the Boltzmann gas. The 
enhancement disappears at higher a s where R/T, and hence the IR cut-off, becomes 
bigger. 

In case (b) and (c) , we introduce an IR cut-off mo by replacing the 1/k 2 -. factor in Eq. 
( TT2l with l/(kj, + TMj)). Thus, tud not only screens the intermediate states but also the 
external states. This is motivated by demanding the optical theorem to be valid in the 
medium, even though it need not be the case when the system exchanges particles from 
a thermal bath. Thus, if the propagator in the loop is screened, then the bremsstrahlung 
gluon is also screened. This very naive treatment gives \kx\ 2 > m 2 D = 0(a s ), which is 
consistent with the first treatment in the a s counting. 



V. NUMERICAL RESULTS AND DISCUSSIONS 



We show in Fig. 2 the comparison between rj computed with gg —>■ g o a 



one (denoted as 

7722) and 7] computed with gg — > gg and gg <-> ggg (denoted as 7722+23) |49(. In computing 
the 23 contribution in case (a) and (b) with the 3-body-like phase space for three gluons, 
we use different treatments of kj< cut-offs: in case (a) we use R = R22 + 2.5i?23 as the 
cut-off, where R23 is self-consistently determined (the blue dashed line in Fig. [1]), while 
in case (b) we use as the regulator. For these two cases we find that adding gg «-> ggg 
reduces 77 by ~ 30% at a s = 10~ 3 where the contribution from gg <-> ggg is about 1/2 
of that from gg <-> gg. The correction is the largest, ~ 75%, at a s = 0.1. This means 
the gg <-> ggg contribution is about 3 times that of gg — > gg. The behavior shown here 
is different from that of XG which shows 7722+23/7722 ~ 1/8 ~ 12.5%, meaning that the 
99 999 contribution is about 7 times as large as gg — > gg, for a wide range of a s 
(a s = 10~ 3 — 0.7). The difference between our result and XG's is largely due to the factor 
2 difference in collisional integrals for the gg <-» ggg process and the BE statistics versus 
the Boltzmann gas approximation used. But we do see the dominance of gg «-> ggg over 
99 ~^ 99 when a s >2x 10~ 3 , as asserted by XG. 

For case (c) with the 2-body-like phase space for three gluons the effect of the 23 
process is about 10-30%, which is close to AMY's result in the whole range of a s . Since 
our result changes dramatically after imposing the soft bremsstrahlung approximation, it 
means this approximation is not as good as previously expected. Thus, it is important to 
go beyond this approximation to obtain an accurate 77. 

In Fig. 3, rj/s as a function of a s is shown for different cases: the LL result t]ll/s of 
Eq. (TIE!) . 7722/s, and 7722+23/s for two different kx cut-offs for the 3-body-like phase space 
and that for the 2-body-like phase space. When a s — >• 0, all these curves should converge 
to the LL result. But at a s = 10~ 3 , we have ln(l/g) = 2.2, which is not large enough to 
dominate the contribution. This is the reason for the deviations of the numerical results 
from the LL one in the current range of a s . However, the agreement between the t/ll/s 
and 7722/s is a good check to our numerical calculations which are carried out by the 
Monte Carlo method for mult i- dimensional integrations. The power of a s dependence of 
these curves are close to (—2) as expected. At a s = 0.1, with both cut-offs for the 3- 
body-like phase space, the full result 7722+23/5 — 1-0 is between 2.7 of AMY 50j and 0.5 of 
XG. At a s = 0.3 and 0.6, we have 7722+23/s — 0.22 and 0.15, respectively, which are larger 




FIG. 2: (color online) 7722+23/7722 shown as a function of a s for the 3-body-like and 2-body- 
like phase space (PS) of three gluons. There are two different treatments of the cut-off of the 
bremsstrahlung gluon momentum fcy for the 3-body-like phase space. 

than 0.13 and 0.076 obtained by XG. It is also interesting to note the good agreement 
using two different cut-offs for the bremsstrahlung gluon momentum. For the 2-body-like 
phase space the correction from the 23 process is small and 7722+23/5 ~ (70% ~ 90%)?722/s, 
which is close to AMY's result. 

In summary, we have calculated the shear viscosity over entropy density rj/s of a gluon 
plasma in kinetic theory. Due to the suppressed contribution to 77 in the gg — > gg forward 
scattering, the gluon bremsstrahlung gg «-> ggg process also contributes at the same order 
(0(a~ 2 )) in perturbative QCD. We find that the gg <-> ggg contribution becomes more 
important to 77 than gg — > gg for q s >2x 10~ 3 for the 3-body-like phase space for the 
three-gluons state. At a s = 0.1, rj/s ~ 1.0 which is between 2.7 obtained by Arnold, 



Moore and Yaffe 



22 



and 0.5 obtained by Xu and Greiner [2l|. Our 77/s is about 2 times 
as large as that of Xu and Greiner for a s > 0.1, largely due to the factor 2 difference in 
collisional integrals for the gg <-» ggg process and the Bose-Einstein statistics versus the 
Boltzmann gas approximation used. We have observed that using as the regulator for 
transverse momentum of the soft bremsstrahlung gluon agrees well with that using the 
rate as the cut-off for the LPM effect in 77 for the current range of a s . In dealing with the 




FIG. 3: (color online) r//s versus a s for (a) the leading- log result in Eq. (|15|) . (b) the result of 
the 22 process only, the full result with 22+23 processes for the 3-body-like phase space (PS) of 
three gluons where the kx cut-off is set by (c) mn or (d) the LPM effect, and (e) the full result 
with 22+23 processes for the 2-body-like phase space (PS). 

23 process it is consistent to implement the soft gluon condition in the energy-momentum 
configuration of the three-gluons state that there is one soft gluon, which results in the 2- 
body-like phase space for the three-gluons state, since we use the Gunion-Bertsch formula 
for the 23 matrix element which is valid only for soft gluon bremsstrahlung. In this 
case we obtain results close to AMY's. To test which is the correct description for the 
phase space of three gluons in the 23 process, or in other words, to test if the Gunion- 
Bertsch formula is still valid for general 3-body-like momentum configurations, a further 
and comprehensive study with the exact matrix element is needed. 
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APPENDIX A: THE CROSS SECTION FOR 23 FROM GUNION-BERTSCH 
FORMULA 



In the center-of-mass frame of 1 and 2, the cross section is written by, 
1 1 /"-q d 3 k 

i=3 



a23 ~ 2^3! 



/ n (2^)3m |Mi2;345|2(27r)V(fci + h ~ h ~ h ~ h) 

= V 2 ^ J d% (c£ + miy m + E *~ E *- ^ S ^ KUT-trY + ^l] 
27 f If q 2 

= J " V ST^ J d2kTdv kU<MT-ly< + miY (A1) 

Since we use the Gunion-Bertsch formula for soft gluon bremsstrahlung, we assume the 
5th gluon is soft, so we made the approximation in the second equality of Eq. (lAip . 

S 4 (h + k 2 -k 3 -k 4 - k 5 ) « 5\h + k 2 -k 3 - h), (A2) 

which means the phase space is dominated by the 22 process. We also used E 3 = E± = 
E 1= E 2 = y/s/2 and 

J d 3 k 3 6(E 1 + E 2 - E 3 - E 4 ) = X - j d 3 q5(E 1 - E 3 ) 

= \J d 2 g T dq z S{E 1 - yJ(E 1 + q z ) 2 + q 2 ) 



d 2 q T ^L= « / d\ T (A3) 
^/E\ -q\ J 

where k3 = ki + q. Note that a factor of 2 is given from the two roots for q z in the 



equation E\ = \J (E~i + q z ) 2 + q\, i.e. q z = —E\ ± \JE\ — q\ which correspond to forward 
and backward solution q z = —y/s, or t = —s, at q? = 0. Eq. ( 1A1I) is 2 times as large 



as that derived in Ref. 
1/2, 



44j. One has to choose the forward scattering and get the factor 



/ d 3 k 3 5(E 1 + E 2 - E 3 - E A ) = ~ / d 2 q T . (A4) 

J forward ^ J 



' forward 

Then the differential cross section from Eq. (lAlj) becomes 



d 2 q T d 2 k T dy 2tt 2 s (q 2 + m 2 D ) 2 k 2 [(q T - k T ) 2 + m 2 ,] 



which reproduces the result in 
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In Fig. 22 of [23j, Arnold has shown the LO ij/s with estimated higher order error from 
two different calculations that should agree at LO but differ at higher order. We see that at 
a s = 0.125, the error is ~20% (4irr]/s = 27 ± 6). The higher order error is more significant 
at larger a s . This error, although significant, is still consistent with an 0(a s ) error. 



